Process tomography of ion trap quantum gates 
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A crucial building block for quantum information processing with trapped ions is a controlled- 
NOT quantum gate. In this paper, two different sequences of laser pulses implementing such a 
gate operation are analyzed using quantum process tomography. Fidelities of up to 92.6(6) % are 
achieved for single gate operations and up to 83.4(8) % for two concatenated gate operations. By 
process tomography we assess the performance of the gates for different experimental realizations 
and demonstrate the advantage of amplitude-shaped laser pulses over simple square pulses. We also 
investigate whether the performance of concatenated gates can be inferred from the analysis of the 
single gates. 

PACS numbers: 03.65.Wj, 03.67.Lx, 32.80.Qk 



Processing information with well-controlled quantum 
systems has the fascinating perspective of being much 
more powerful than classical computers for certain ap- 
plications. A promising candidate for the experimen- 
tal realization of quantum computing are strings of ions 
stored in linear Paul traps [l| as recently demonstrated 
by various key experiments, including the preparation of 
multi-particle entangled states 0, 0] , quantum teleporta- 
tion 0, and quantum error correction . Quantum 
information processing depends on the ability to imple- 
ment single qubit rotations and most importantly an en- 
tangling two-qubit quantum gate 0, El E3 ■ Proper 
characterization and understanding of the action of gate 
operations and their imperfections is of vital importance 
in order to successfully apply them in complex computa- 
tions. 

Generally, the implementations of quantum gates are 
imperfect due to decoherence and various systematic er- 
ror sources present in experimental setups. A proper 
description of such an operation which accounts for the 
possibly non-unitary evolution of the qubits is provided 
by quantum process tomography [TlLll2j . Process tomog- 
raphy has already been applied for characterizing quan- 
tum gates in NMR and linear-optics quantum comput- 
ing Here, we show that process tomography 
is a valuable tool for comparing different ion trap quan- 
tum gate implementations and optimizing the experimen- 
tal parameters. This way, we were able to improve our 
controlled-NOT (CNOT) gate fidelity from 71 % to 
almost 93 %. Moreover, the action of two successively 
applied gate operations is investigated and compared to 
the predictions from the single gate tomography result. 

We realize entangling gates between 40 Ca + ions held in 
a linear trap |l6|. Quantum information is stored in su- 
perpositions of the \S) = 5i/ 2 (to = —1/2) ground state 
and the metastable \D) = D 5 / 2 ( m = —1/2) state and is 
manipulated by laser pulses at a wavelength of 729 nm 



exciting the electric quadrupole transition between those 
states. A focus size smaller than the inter-ion distance 
and precise control of the focus position allows us to ad- 
dress single qubits. Detection of the qubit's quantum 
state is achieved by scattering light on the S±/2 <-> P1/2 
dipole transition and detecting the presence or absence 
of resonance fluorescence of the individual ions with a 
CCD-camera. The oscillations of the ions in the har- 
monic trap potential are described by normal modes and 
give rise to sidebands in the spectrum of the S1/2 <-> -D5/2 
transition. For coherent state manipulation, only the 
quantum states \n) of the axial center-of-mass mode at 
a frequency w z = 2ir ■ 1.36 MHz are relevant. Here, n 
denotes the number of vibrational quanta. Quantum in- 
formation processing is implemented by (a) laser pulses 
on the carrier transitions \S,n) <-> \D,n) realizing single 
qubit operations on the ion qubits and (b) laser pulses 
on the first blue sideband inducing transitions between 
the states \S, n) and \D,n + l) which connect the internal 
state of the ions and the state of the vibrational mode. 
The latter operation allows us to implement an entan- 
gling interaction between ion qubits |]| . A more detailed 
account of our experimental setup can be found in Ref. 

M 

Cirac-Zoller controlled-NOT gate operations between 
two ion qubits are implemented by the pulse sequences 
shown in Tab. [I] In both sequences, the quantum state 
of the control qubit 2 is first mapped to the vibrational 
mode (SWAP I), which is cooled to its ground state prior 
to the operation. Then a CNOT operation is performed 
between the vibrational mode and the target ion qubit 1. 
Finally the state of the vibrational mode is mapped back 
onto the control qubit (SWAP II), restoring its quantum 
state and returning the vibrational mode to its ground 
state. Both pulse sequences differ only in the way the 
phase gate sequence is implemented. The ideal unitary 
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evolution realized by the first pulse sequence (A) is Q 
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where the matrix is written in the product basis 
{\DD), \DS), \SD\, \SS}} and expressed in terms of the 
Pauli operators X, Y, Z and the identity I. In this 
sequence the state of target ion 1 is flipped whenever 
control ion 2 is in state \D) (the order of the ions is 
|ion 2, ion 1)). The unitary evolution of the second pulse 
sequence (B) is 
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Here the state of the target ion 1 is flipped whenever 
the control ion 2 is in \S). This pulse sequence shows 
the desired unitary evolution if an additional rotation 
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applied. If this CNOT pulse sequence is embedded in 
a larger algorithm, the additional rotation can be taken 
into account by shifting the phase of every subsequent 
pulse by A0 = —1/^/2 ■ it on the control ion and by 
• 7r on the target ion. As can be seen from 
Tab. [Q pulse sequence (B) is shorter in terms of total 
length of the sideband pulses than sequence (A). 



Description 


Seq. (A) 


Seq. (B) 


SWAP I 


Jtf(ir,0) 


Jtf(ir,0) 


Ramsey I 


-Ri(tt/2,0) 


Ri(tv/2,0) 


Phase gate 


R+^,n/2) 
R+(tv/V2,0) 
Bf 

R+(n/V2,0) 


R+ (n/2,n) 
Rf(V2n, tt/2) 
flf(ir/2,0) 


Ramsey II 


Ri(ir/2,n) 


fli(jr/2,(l/>/2-l)w) 


SWAP II 


i?+(7T,7r) 


J?+(7r,7r) 



TABLE I: Two sequences of laser pulses for implementing a 
CNOT gate operation. Laser pulses applied to the i-th ion 
on the carrier transition are denoted by Ri(8,<f>) and pulses 
on the blue sideband transition by Rf(6, 4>), where 6 = Q • t 
denotes the pulse area in terms of the Rabi frequency Q, the 
pulse length t and its phase <f> [Ts|| . 

Due to systematic imperfections and decoherence the 
actual evolution in our experiment will deviate from the 
ideal unitary evolution given in (JJJ and (J2J- An impor- 
tant systematic error in our setup is imperfect address- 
ing. Every laser pulse which addresses one of the ions 



also slightly affects the neighboring ion qubits, due to 
the finite size of the laser beam focus. In terms of the ra- 
tio of Rabi frequencies between the addressed and the 
neighboring ions, this error is on the order of 2-3 %. 
Furthermore, due to decoherence the output state after 
application of an operation will in general be a mixed 
state. The major sources of decoherence in our experi- 
mental setup are fluctuations of the laser frequency and 
the ambient magnetic field [lflj . 

The experimentally realized quantum gate including 
error sources is properly described by a completely posi- 
tive map £. For an input state p, the output state £{p) 
can be written in the operator sum representation as 
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where N is the number of qubits and the A m are operators 
forming a basis in the space of 2^ x 2 Ar -matrices. The 
process matrix x contains complete information about 
the investigated process including the influence of the 
environment on the qubits. The matrix \ can be ex- 
perimentally obtained by employing quantum process to- 
mography. This procedure requires 4^ input states pi n> \ 
which are linearly independent. For every input state, 
the output state £(pi n ,i) has to be determined by quan- 
tum state tomography. From this set of data the process 
matrix \ can be obtained by inverting the relation in y| . 
However, to avoid unphysical results caused by quantum 
noise in the measurement process, we employ an iterative 
maximum likelihood algorithm |2(| in order to find the 
physical process £ which most likely generated the mea- 
sured data set. We choose the products of the single qubit 
states 1^) = \S), |^ 2 > - \D), |^ 3 ) - (\D) +i\S))/y/2 and 
1^4) = (\D) + \S))/V2&sthe 16 input states necessary for 
a tomography of our two qubit quantum gates. Quantum 
state tomography of a two qubit system then requires 
measurements in nine different product state bases |l8| . 
This results in a total of 16 x 9 = 144 different measure- 
ment settings. Every expectation value is determined 
through 100-250 individual experiments at a rate of 50 
experiments/s which requires 5-12 minutes of measure- 
ment time to gather all data necessary for the estimation 
of x- 

Quantum process tomography was carried out for the 
operations implemented by pulse sequences (A) and (B) 
shown in Tab. [I] The duration of the complete CNOT 
pulse sequence was Toate = 615 ps for sequence (A) and 
Tbatc = 502 ps for (B), which is mainly determined by 
the Rabi frequency on the blue sideband at Obsb = 2n ■ 
4.4 kHz. The resulting estimated process matrix Xcnot 
for pulse sequence (B) is shown in Fig. ^ 

The process matrix allows us to calculate various mea- 
sures which characterize the performance of the gate op- 
eration. We can directly calculate the process fidelity 
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FIG. 1: Process matrix obtained by process tomography of a gate operation implemented by pulse sequence (B) in Tab.[TJ The 
absolute value, real part and imaginary part of xcnot are shown in a), b) and c) respectively. The matrix is expressed in terms 
of the products of the identity / and the Pauli operators X, Y and Z. In order to compensate for the rotation Uz the phase 
of all tomography pulses was shifted appropriately. 



F p = tr (xid ■ Xcnot), which is the overlap of our exper- 
imentally obtained xcnot with the ideal process matrix 
Xid derived from the unitary evolution of the gate given 
in (JJJ and (J2J), respectively. Furthermore, using Eq. (0 
we can predict the output state of the experimental gate 
operation for an arbitrary input state. This enables us 
to investigate the gate performance for a large number of 
numerically generated input states, similar to the analy- 
sis done for an optical CNOT gate in Ref. [lif. We do 
this by analyzing the calculated output states in terms 
of their overlap F — ( , 0id,out|£(/o)|V ; id,out) with the ideal 
output states (f/'id.out), their normalized linear entropy 
S\in = 4/3-tr (l — £(p) 2 ) and the change in entanglement 
from the input to the output states given by the change in 
the concurrence, i.e. the difference AC = C(£(p)) — C(p) 
|2l) . For 5 • 10 4 randomly chosen pure input states drawn 





Soq. 


F p in % 


Fmean in % 


Sim 


max(AC) 


Description 


1 


A 


88.8(7) 


91.0(6) 


0.20(1) 


0.86(2) 


single gate £ C not(A) 


2 


B 


90.8(6) 


92.6(6) 


0.17(1) 


0.84(2) 


single gate £ C not(b) 


3 


A 


87.7(7) 


90.1(6) 


0.21(1) 


0.81(2) 


with pulse shaping 


4 


A 


75(1) 


80(1) 


0.39(2) 


0.70(3) 


no pulse shaping 


5 


AA 


79(1) 


83.4(8) 


0.34(2) 




£2xCNOT(A) 


6 


A 


82.8 


86.2 


0.30 




£cNOT(A) £cNOT(A) 


7 


BB 


72(1) 


77.4(8) 


0.41(1) 




^2xCNOT(B) 


8 


B 


79.8 


83.8 


0.34 




^CNOTfB) £cNOT(B) 



TABLE II: Performance of gate operations for different ex- 
perimental realizations. The mean fidelity F meall , the mean 
linear entropy Sn n and the maximum change in entanglement 
max(AC) were inferred from an ensemble of 5 ■ 10 4 random 
states. The quoted errors are due to quantum projection noise 
|18|. For the results in rows 1,2,5-8 the blue sideband Rabi 
frequency was set to flsss = 2-7T-4.4 kHz. In rows 3,4 a higher 
Rabi frequency of P2bsb = 2-7T • 5.3 kHz was chosen. Rows 5,7 
contain the results for the tomographies of two concatenated 
gates. Additionally, rows 6,7 show the results predicted from 
the single gate analysis. 



from the Haar measure on the unitary group U(A), we 
carry out such an analysis using the results of the process 
tomography for the two types of quantum gates described 
above. We characterize the gate performance by calculat- 
ing the mean fidelity i^mean, the mean linear entropy S\\ n 
and by searching for the maximum increase in entangle- 
ment max(AC). The results for a single gate operation 
are shown in Tab. ITU (rows 1,2). Pulse sequence (B) shows 
a slightly better performance than sequence (A) probably 
due to the 20 % shorter phase gate sequence, which re- 
duces the influence of decoherence. Furthermore, errors 
due to imperfect addressing will partially cancel in phase 
gate sequence (B), since the two 7r/2-pulses are applied 
with a relative phase of 7t. 

We used this kind of analysis to assess the influence of 
amplitude pulse shaping on the gate performance. Excit- 
ing the blue sideband transition with square laser pulses 
causes off-resonant excitation of the carrier transition de- 
grading the performance of the quantum gate. If higher 
gate speeds were to be achieved by using higher side- 
band Rabi frequencies [2^, this effect would become in- 
creasingly harmful. However, we largely suppress off- 
resonant excitations by adiabatically switching on and 
off the laser pulses 23]. We demonstrated this by first 
carrying out a process tomography of gate pulse sequence 
(B) for Tg a t c = 520 /is using shaped pulses, that were adi- 
abatically switched on and off with a rise and fall time 
of 5 ps (the standard setting for all reported results). 
Then we carried out another gate tomography using sim- 
ple square pulses. As can be seen from the results in Tab. 
ITTI (rows 3,4), the use of shaped pulses considerably im- 
proves the gate performance compared with the result 
with square pulses. 

Quantum algorithms will generally contain multiple 
quantum gates, which are successively applied to a qubit 
register. The question arises whether the performance of 
such a series of quantum gates can be inferred from the 
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knowledge of the single-gate performances in experimen- 
tal implementations of an algorithm. We investigated 
this issue by comparing the result of a process tomog- 
raphy of two concatenated CNOT gates with the pre- 
dictions inferred from the single gate tomography, for 
the same set of experimental parameters. Ideally, two 
concatenated CNOT yield the identity, (J7cnot) 2 = I- 
The measured process matrix X2xCnot is shown in Fig. 
El As expected, the dominant element of X2xCnot is the 
77, //-element. We determined the process fidelity, mean 
fidelity and mean entropy of the operation described by 
X2xCnot- For comparison, we calculated the same quan- 
tities for the process £cnot ° £cnot, where £cnot was 
obtained from the single gate process tomography. The 
results are shown in Tab. [H] (rows 5-8). In general, 
one would expect both methods to yield the same re- 
sults if the dynamics of interaction between the qubits 
and the environment was Markovian, thus producing un- 
correlated errors in both gates. In our experiment, we 
attribute the observed discrepancy to low- frequency mag- 
netic noise giving rise to magnetic fields that are constant 
over the course of the double gate sequence but vary from 
experiment to experiment. In addition to characterizing 
the performance of gates within larger blocks, concate- 
nating quantum gates might be useful for amplifying tiny 
errors in high-quality gates to a measurable size. 




FIG. 2: Absolute value of the measured process matrix result- 
ing from two gate operations successively applied to a pair of 
ion qubits. Ideally the process matrix should only contain 
the //, //-element. The height of the //, //-element in the 
measured process matrix of 79.3 % directly gives the process 
fidelity. 

The presented work demonstrates quantum process to- 
mography to be a valuable tool for assessing the perfor- 
mance of the fundamental operations of ion trap quan- 
tum computers. Since a complete tomographic data set 
can be taken in a comparably short amount of time, it is 
easily possible to compare the performance of quantum 



gates for different pulse sequences and experimental pa- 
rameters. This has been illustrated by showing that adia- 
batically switching on and off laser pulses instead of using 
rectangular shaped pulses make higher gate speeds pos- 
sible while preserving high gate fidelities. This technique 
and careful optimization of the experimental parameters 
using process tomography helped to significantly improve 
the gate fidelity compared with earlier results jjJLLil The 
results of the concatenated gate tomography demonstrate 
the importance of analyzing experimental quantum gate 
implementations not only as isolated objects but also 
within a larger gate sequence. We expect that this tech- 
nique will have considerable impact on estimating the 
overall performance of future quantum computers. 
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